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SUMMARY 


A new method is developed for solving the lifting-surface equation for 
thin wings. The solution requires the downwash equation to be in the form of 
Cauchy integrals which can be interpreted as a vortex lattice with the posi- 
tions of the vortices and control points dictated by the finite sum used to ap- 
proximate the integrals involved, ban's continuous loading method is employed 
for the chordwise integral since it properly accounts for the leading-edge 
singularity, Cauchy singularity, and Kutta condition. Unlike Lan, the spanwise 
loading is also continuous and the Cauchy singularity in the spanwise integral 
is also properly accounted for by using the midpoint trapezoidal rule and the 
theory of Chebychev polynomials. This tecnuiq^’^ yields the exact classical so- 
lution to Prandtl’s lifting-line equation. The solution to the lifting-surface 
equation for rectangular wings was found to compare well with other continuous 
loading methods, but with much smaller computational times, and it converges 
faster than other vortex lattice methods. 


INTRODUCTION 


The vortex lattice method has proven to be a useful technique for calcula- 
ting the aerodynamic characteristics of complete configurations as well as 
wings. In the conventional vortex lattice method (VLM), the planform is di- 
vided into a number of elemental panels, and a horseshoe vortex is placed at 
the local quarter-chord of each panel. The boundary condition is satisfied at 
the local three-quarter chord of each elemental panel (called control points) 
by requiring the flow to be tangent to the surface there. The strengths of the 
horseshoe vortices are determined by solving the matrix equation formed from 
the tangent-flow boundary conditions. Then, the aerodynamic characteristics 
are calculated by summing the results from each elemental panel. A complete 
description of the conventional vortex lattice method is given by Margason and 
Lamar in reference 1. 

Although reasonable results are obtained by the conventional VLM, Lan (ref. 
2) listed the following deficiencies: 1) The method used to compute the in- 

duced drag implies that the leading-edge thrust is distributed over the chord 
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instead of being concentrated at the leading edge. 2) The predicted pressure 
distribution is not accurate near the leading edge. 3) The convergence of so- 
lutions is slow with respect to the number of panels used. In addition, the 
Kutta condition is not explicitly satisfied. Hough (ref. 3) found some im- 
provement by using a 1/4 lattice width inset at the wing tips. 

Inn (ref. 2) developed an ingenious method for thin, two-dimensional air- 
foils by using the midpoint trapezoidal rule and the theory of Chebychev poly- 
nomials to reduce the downwash integral to a finite sum. This method gives the 
exact lift, pitching moment, leading-edge suction, and pressure difference at 
a finite number of points; and the Kutta condition is satisfied at the trailing 
edge. A more detailed description is given below. Ian also developed a quasi- 
vortex lattice method for finite wings by using his two-dimensional method for 
the continuous chordwise vortex distribution but a stepwise constant vortex 
distribution in the spanwise direction. The results showed an improvement over 
those calculated by the conventional vortex lattice method. 

This paper develops a new vortex lattice method which uses Lan*s continu- 
ous chordwise vortex distribution but, unlike Lan, a continuous spanwise vortex 
distribution also. Although the vortex distributions are continuous, the 
method is easily interpreted as a vortex lattice method in which the arrange- 
ment of horseshoe vortices and control points are determined from the finite 
sum used to approximate the downwash integral of lifting-surface theory. In 
order to understand the development of the present method, ban's two-dimensional 
theory is reviewed first, and then the present method is applied to Prandtl's 
lifting-line theory before developing the method for lifting-surface theory. 


SYMBOLS 


A aspect ratio 

b wing span 

c wing chord 

c^ sectional lift coefficient 

c sectional moment coefficient about leading edge 

m 

sectional leading-edge thrust coefficient 

C far-field induced drag coefficient 

°i 

C near-field induced drag coefficient 

■^ii 

wing lift coefficient 
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lift-curve slope, per radian except when noted otherwi^ie 


w: g pitching moment coefficient about leading edge 


pitching-moment curve slope, per radian 


leading-edge suction parameter 


wing leading-edge thrust coefficient 


far-fleld spanwlse efficiency factor, C. /C TT A 


near-field efficiency factor, C /C tt A 

L D.^ 


parameter defined by eq. (26) 


number of trailing vortices over whole wing span 
number of spanwise control points over whole span 
summational Integer 

number of chordwise vortices and control points 
National Aerospace Laboratory, Netherlands 
wing planform area 
vortex lattice method 
freestream velocity 

downwash velocity, referred to and positive upwards 
chordwise coordinate measured from leading edge in direction of 
sectional and wing aerodynamic center locations, respectively 
spanwise coordinate, positive to the right 


vertical coordinate of mean camber line 
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angle of attack 

nondimensional circulation per unit chord 
circulation 

difference between lower and upper pressure coefficients (AC = 

P 

transformed chordwise coordinate, see eq. (2) 
transformed spanwise coordinate, see eq. (12) 


Subscripts : 


chordwise control point, see eq. (7) 
spanwise control point, see eq. (15) 
chordwise vortex position, see eq. (6) 
spanwise trailing vortex position, see eq. (14) 
evaluated at spanwise position - pTf/M 


LAN’S TWO-DIMENSIONAL THEORY 


For thin airfoils, the downwash equation is 


1 f '**1 


27T J 

o 


The integral on the right side is of the Cauchy type. Transform the x coordi- 
nate by 

x/c = (1 - cos 0)/2 (2) 

and use the following result from airfoil theory (ref. 4) 




J cos 0^ - 


cos 0 


to write eq. (1) as 
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T 


W. 


1 I sin d0^ ^ [y(0j^) sin 0^^ - Y(0) sin 0]d0j^ 

2tt j cos 0j^ - cos 0 27T 1 cos 0, - cos 0 


cos 0j^ - cos 0 


(4) 


N 

1 


, 4, cos 0, - cos 0, 
k=l k i 

when the vortex positions are 

A - (2k - 1)TT 
®k 2N 

and the control points are located at 

iTT 


polynomials to show that 


-n2 

for 

i = 0 


0 

for 

i 0, N 

(5) 


for 

i = N 



} k“lj 


®i “ T ’ i “ 0, 1 N 


( 6 ) 


(7) 


Then the integral in eq. (4) can be reduced to a finite sum by using the mid- 
point trapezoidal rule (ref. 5) and eq. (5) to obtain 

, ^ N Yi sin 0, - y. sin 0, 

1 7T r 'k k i i 

N , ^ 


'^i 2ir N , cos 0, - cos 0. 

k*l k i 


N 


2N , cos 9, - cos 0, 
k=l k i 


= _ i y 

OM L 


- y lim y( 0) sin 0 , i = 0 

^ 0-H) 




, i 0, N (8) 


^ lim Y(0) sin 0 , i = N 

^ 0-^ 


However , 


lim Y(0) sin 0 = 4 
0-K) 


(9) 


where C„ is the leading-edge suction parameter and since the Kutta condition re- 
quires that Y(Tf) = 0, 


lim Y(0) sin 0=0 (10) 

0-»TT 


Unlike the conventional vortex lattice method, the Cauchy singularity, 
leading-edge square-root singularity and the Kutta condition are properly ac- 
counted for in this method. Equation (8) can be solved with 1 0 to obtain 



t 


\ 
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the N values of Yu» then the leading-edge suction parameter can be computed 
by using eq. (8) with 1=0 (control point at the leading edge). Figure 1 il- 
lustrates the positions of the vortices and control points by the "semicircle 
method" for N = 2 . With only one vortex (N = 1) , the exact lift and leading- 
edge suction are obtained, and the Kutta condition is satisfied. With two or 
more vortices the exact pitching moment is obtained in addition to the above 
properties and the calculated values of are exact. It can be shown that the 
remarkable accuracy of this method is due to eq. (5), which is similar to the 
integral result given by eq. (3) and used in exact thin airfoil theory. 


PRANDTL'S LIFTING-LINE THEORY 


Before attacking the lifting-surface equation, the present method will be 
developed for Prandtl's lifting-line equation to compare the spanwlse lift dis- 
tribution with the classical solution (ref. 4). The lifting-line equation is 
given by ref. 6, 


r = 7T V c a - 


dr 

(y - y^) 


This equation also has a Cauchy integral on the right side, and thus Lan's 
airfoil technique is applicable. Transform the spanwise coordinate by 

y = - y cos (f> (1 

and replace the downwash integral in eq. (11) with the midpoint trapezoidal- 
rule summation to get 


4tt V 


b/2 

1 f dr 


X_TT2 r 

T V M h I 


(dr/d4>^) 


dy^ (y^ - y^) 4tt M b cos - cos 


This equation represents the downwash due to M trailing vortices of strength 
-(dr/d(()j^)^ (tt/M) located at 

_ (21 - pTT 0 _ 1 w f-, 

" 2M , X — 1, . . . , M (J 


with control points located at 


, M-1 
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Now the conventional vortex lattice arrangement Is shown in figure 2*, whereas 
the present arrangement is shown in figure 3. Note in particular that the tip 
vortices extend to the wing tip in the conventional method, whereas eq. (lA) 
determines them to be inset (see fig. 3), which agrees with Hough's results 
(ref. 3). 

Since (dr/diJij^)j, is needed in eq. (13), it is convenient to represent the 
circulation by Multnopp's interpolation formula 


M-1 M-1 


“ M I I 

P=1 P n=l P 


(16) 


where 




= £2L 


p “ 1, ..., M-1 


p M 

Equation (16) io based on the following orthogonality property (ref. 7) 


(17) 


^ M-1 I ^ 

[ sin n6 sin n(j) = < 
n=l ^ 1 0 


for p = j 


for p j 


(18) 


Substitute eq. (16) into eq. (13) and then the downwash becomes 

M M-1 M-1 sin nij)^ n cos nc))^ 


w. - 


I I r„ Z 


J V 1-1 p-l P n-l *1 ■ ♦j 

00 *■ j 


(19) 


However, this result can be simplified by using eq. (5) to derive the important 
summation below 


t 


( ‘ 
( 


H 

M 


M 

I 


£=1 


cos n(()^ 

cos (f)^ - cos (p 


j 


TT sin nd). 
^ 


sin (j) 


J 


( 20 ) 


This equation is similar to the following Integral result used in thin airfoil 
and lifting-line analyses (ref. A) 


J 


cos nc|)j^ d(J)j^ 
cos (j)j^ - cos <(i 


TT sin n(j) 
sin (p 


( 21 ) 


When eq. (20) is used in eq. (19), the downwash reduces to 


In some conventional techniques the spanwise length of each vortex is uniform, 
but it is reported to have little effect on the results. 
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1 

V M b 


M -1 M -1 sin nd) n sin n({). 

y r y -E ^ 

p=l P n=l 


Finally, substitute eq. (22) for the downwash in eq. (11) and apply the result- 
ing equation at the spanwise locations (j)j(j = 1, M-1) to obtain a matrix 

equation for the (M-1) values of Fp. Then, the lift, pitching moment, and in- 
duced drag can be calculated by reducing the spanwise integrals to a finite sum 
through the midpoint trapezoidal rule. However, the midpoint trapezoidal rule 
gives exactly the same result for the spanwise integrals as the integral itself 
when eq. (16) is used for F. To prove this assertion, consider the spanwise 
integral for the lift. Using eq. (16), the exact integral result is 


Ct V ^ 
L 00 b 


11— X 

“ T7 J ^ sin (f) 

M ^ n 


M p 

p=l 


whereas the midpoint trapezoidal rule gives 


y" M M-1 M-1 M-1 

“ Fdy=^— y 1 sin ncf) sin n((). sin ~ ^ 1 ^ ^ 

b I n T 1 P 1 P ^ ^ M ^1 P t 


which is the same as the exact integral result above. In obtaining this last 
result, the following orthogonality property was applied 


J sin n(p^ sin (p^ = 


T for n = 1 


for n ^ 1 


The integrals for the induced drag are handled in a similar fashion. 


The remarkable feature of the present method is that the results are iden- 
tical to the classical solution of Prandtl*s lifting-line equation when a 
finite number cf terms is used in the Fourier series for F (ref. 4). The suc- 
cess of the present method is attributed to the location of the spanwise vor- 
tices and control points, the summational result of eq. (20), and the accuracy 
of the midpoint trapezoidal rule for the spanwise integrals. 

Figure 4 Illustrates the convergence of this method compared with the con- 
ventional VLM for Cx^ of rectangular and elliptical planforms with an aspect 
ratio of 2 tt. Prandtl’s lifting-line theory requires trailing vortices and 
spanwise control points, but no chordwise control points are needed because the 
theory assumes the downwash is constant in the chordwise direction. Therefore, 
the control points are placed on wue **bound” vortex for both methods. In this 
way, the accuracy of the spanwise vortex arrangement can be tested without the 
influence of the location of chordwise control points. Figure 4 shows that the 
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conventional VIM converges slowly and Ci^ does not appear to approach the cor- 
rect limit when the curve Is extrapolated to an Infinite number of trailing 
vortices, M^(l/M^). On the other hand, the present method converges very 
quickly and approaches the correct limit. As In the classical solutions, only 
one horseshoe vortex (M»2) is needed to obtain the exact for the elliptical 
planfonn in Prandtl's llftlng.llne theory. Figure 5 illus^ates the effect of 
the number of trailing vortices on the spanwlse efficiency factor used In the 
Induced drag. Here again, the conventional VIM converges slowly and does not 
appear to approach the correct limit, whereas the present method converges 
quickly and approaches the correct limit. 


LIFTING-SURFACE THEORY 


For simplicity, the present method is developed here for rectangular wings. 
The downwash equation from lifting-surface theory is usually given by one of 
the two following integrals (ref. 6) 




Y(Xj, y^^) 

s 'S’ - >'l>^ 


1 + 


(X - Xj^) 


Ath - Xj^)^ + (y - 


dxj^ dyj^ (23) 


or 


w(x,y) 


-L f f At 1 
i J 9yi (y - Yi) 
s 


A 


X - 


1 + 


Xj^)^ + 


(y - y^)' 


(x - Xj^) 


dXj^ dyj^ 

(24) 


Equation (23) contains the Mangier-type integral, and therefore, is not suitable 
for the present method. Equation (24), however, contains Cauchy-type integrals 
and is therefore in the form to apply a combination of Lan's two-dimensional 
method for the chordwise integration and the lifting-line method developed above 
for the spanwlse integration. Note that the integrand of eq. (24) represents 
the downwash at (x,y) due to half of a horseshoe vortex as shown below. 




I 


With y * -b cos 0/2 and x/c ■ (1 - cos 0)/2, replace both Integrals In eq. (24) 
with the midpoint trapezoidal-rule sum to get / 


^ IL IL ~ y y 9 y i 

''i.J “ " An M N b (cos 




<|)^ - COS (|)j) 


-2NC„ , i - 0 


0 . 11*0 


where 


''ijkl! - 1 * 


i 2 2 2 

Hccs 0^ - cos 0^) + A (cos - cos 0.) 


cos 0, - cos 0. 
k i 


The control points are located at 


0^ = ^ y i = 1, N (chordwise) 


0j="^,j=l, M-1 (spanwise) 


and the “joint** in the horseshoe vortices are located at 

> k = 1, N (chordwise) 


_ (2& - 1)1T 


, il = 1, M (spanwise) 


The positions of the horseshoe vortices and the control points are Illustrated 
in figure 6 by the **semicircle method** for two chordwise vortices (N=2) and 
four trailing vortices (M“4). As in eq. (16), represent the spanwise variation 
of Yk(^)» chordwise position 0y^, by Multhopp’s interpolation formula, 

M-1 M-1 

" M ^ k ^ 

p=l n=l ^ 

where Yp^k ‘Represents the unknown circulations per unit chord at and 0p * pw/M. 
Substitute eq. (31) into eq. (25) to obtain the final f jnt of the downwash as 


= - ^ n Sin iK^p cos sin 9^^ 

2bn\ £>1 k-1 p=l n=l h - 


-2NCg , i = 0 


0 ,1^0 
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The tangent-flow boiindary condition for thin wings requires that 

rhz \ 


l.j 


where z^(x,y) is the shape of the mean camber line. The N(M-l) values of 
are calculated by solving the matrix equation formed by applying f.q. (32) for 
1 0 at the chordwlse and spanwlse control points given by eqs. (27) and (28). 

•* Then after Yp,k calculated, the (M-1) leading-edge suction parameters 

Cgj can be computed by successively applying eq, (32) with 1 ■ 0 (control point 
ac the leading edge) at the spanwlse positions j > 1, M-1. Regardless of 

the number (N) of chordwlse vortices used, there Is always a control point at 
the trailing edge which satisfies the Kutta condition, and another control point 
at the leading edge which gives the leading-edge suction parameter. If desired. 

The sectional and wing aerodynamic characteristics may now be calculated 
by using the midpoint trapezoidal rule to reduce the Integrals to finite sums, 
as shown below (for rectangular wings) 


2T 

c 

2 f 


- J J YpCxj) djCj = 


0 

b/2 

/ tr M-1 

“ 1 

H 'V® ”2ff \ 

-b/2 

' p=l 


(Cj)„ sin (p 
X, p p 


2 1 ir ^ 

■ ■ J J ’‘l ■'’‘l “ • 2N Jj ’^p.k <1 - V 


' M-1 

Cm - I dy/sc . i j 

, ' P=1 


-b/2 


(c„) sin * 
m p p 


C‘ M-1 M-1 

^ I ^ I ^ 

”*1.1 pii " I 


M-1 1 ' 

I K 



(Al) / 


<s>j ■ 2” 



The spanwlse loading can be made coni.i'’uous by eq. (31), and the chordwise 
loading can also be made continuous by fitting Cgj and Yj chordwise 

loading functions for thin airfoil theory (ref. * 


RESULTS FOR RECTANGULAR WINGS 


For one horseshoe vortex and one control point (N*l, M“2) , the present 
method yields 



7T A 


1 + /l + A^/7 


and 



C 


2 

L 


ttA 


These results give the correct limit as \-+0, but just as Lan found for airfoils, 
at least two chordwise vortices are needed to get an accurate pitching moment. 

Table 1 gives a detailed comparison cf the results of the present method 
with those of several other methods for a flat A * 2 rectangular wing. The 
methods chosen for comparison are the continuous loading method of the National 
Aerospace Laboratory of the Netherlands (NLR) presented in ref. 8, Lan's quasi- 
vortex lattice method (ref. 2), the conventional vortex lattice method of 
Margason and Lamar (ref. 1), and Wagner's continuous loading method (see ref. 2). 
In the NLR method (M-1) spanwise loading functions are applied but 8M spanwlse 
integration points are used. Therefore, :he results from this method are used 
as a base for comparison purposes. Table 1 shows that the present method yields 
more accurate overall aerodynamic characteristics than either the conventional 
VLM or Lan's quasi-vortex lattices. This table also shows that the spanwlse 
variation of the sectional lift coefficient compares within four significant 
figures to those of the NLR method. The spanwise variation of the sectional 
aerodynamic center also compares well except near the wing tip. When the nunt>er 
of chordwise vortices was increased from N ■ A to N ■ 6 in the present method, 
these differences decreased considerably. Figure 7 compares the present method 
with the NLR method for the chordwise loading at midspan on this same A ■ 2 
wing. Again, the results compare quite well. 

The computational time required for the results in Table 1 was 22 minutes 
for the NLR method on a CDC 3300 computer; Lan's method required one aiinute on 
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the Honeywell 635 computer} Wagner's method used about three minutes; and 
the present method required less than ten seconds on an IBM 37G/165. Therefore, 
the present method is as economical as the VLM for the same number of vortices, 
but it Is generally more economical when one considers that a smaller number of 
vortices can be used to achieve the same accuracy as the VLM- 

The effect of the vortex lattice arrangement on the convergence of the 
lift-curve slope is presented In figure 8 and compared with the conventional 
VIM for flat rectangular wings with A ■ 2, 4.5, and 7. Tills figure Illustrates 
again the slow convergence of the conventional VLM and the rapid convergence of 
the present method. For all three wings, the present method gives good ac- 
curacy with only two chordwlse vortices, but more spanwise vortices are needed 
for the A ■ 7 wing (M/2 « 10) than the A ■ 4.5 wing (M/2 = 5) or A ■ 2 wing 
(M/2 “ 2). Results for the pitching moment, aerodynamic center location, and 
Induced drag were found to converge even faster than the lift-curve slope, 
therefore they are not presented. 

The Prandtl-Glauert rule can be easily applied to the present method to 
Include subsonic compresslMlity effects. 


APPLICATIONS TO OTHER CONFIGURATIONS 


Flaps and ailerons may be added to the wing by using an approach somewhat 
similar to that of Lan (ref. 2). For the chordwlse integration, the Interval 
from the leading edge of the wing to the flap leading edge is transformed into 
[0,7rl by the "semicircle method", and the interval from the flap leading edge 
to the trailing edge Is also mapped into t0,ir] by another "semicircle". The 
same technique can also be applied in the spanwise direction. 

Application of the present method to tapered and/or swept wings requires 
additional considerations. Care must be exercised so that the chordwlse vortex 
and control points at one spanwise location match those at another spanwise 
position in order to evaluate the Cauchy Integral properly. These configura- 
tions are presently being studied along with non-planar wings. 


CONCLUDING REMARKS 


A new method is developed for solving the I' f ting-surf ace equation for 
thin, subsonic wings. The downwash equation is written as Cauchy-type inte- 
grals for the chordwlse and spanwise directions. They can be interpreted as a 
lattice of horseshoe vortices and the positions of the vortices and control 
points are determined by the finite sum used to approximate the integrals. 

Lan's two-dimensional method is used for the chordwlse integral since it pro- 
perly accounts for the leading-edge singularity, Cauchy singularity, and the 
Kutta condition. For the spanwise integral, Multhopp's interpolation formula is 
used in conjunction with the midpoint trapezoidal rule and the theory of 
Chebychev polynomials. This method properly accounts for the Cauchy singularity 
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and yields the classical solution to Prandtl’s lifting-line equation. The nu- 
merical method for evaluating the chordwioe and spanwise integrals is much 
simpler and quicker than other continuous loading methods. 

The chordwlse and spanwise methods are combined ^ o obtain a continuous 
loading solution to the lifting-surface equation. The algorithm for the rectan- 
gular wing gives results which compare well with other rcatinuous loading 
methods, but with much smaller computational times. In addition, it converges 
faster and is more accurate than other vortex lattice methods. 

For rectangular wings, the vortex lattice arrangement dictated by the pre- 
sent method differs from the conventional VLM in that the chordwise positions 
of the vortices and control points do not follow the usual 1/4 -3/4 rule. 

There is always a control point at the trailing edge, which allows the Kutta 
condition to be satisfied, and a control point at the leading edge which yields 
the leading-edge suction parameter. The spanwise vortices determined by the 
present method are not uniformly spaced, and the tip vortices are inset from 
the actual wing tip. This vortex lattice arrangement gives better results than 
other VLM^s for rectangular wings. Other wing planforms require additional 
considerations, and they are presently being investigated. 
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TABLE 1. RESULTS FOR RECTANGULAR PLANFORM 

A = 2 



Present 
N=4, M=16 

NLR 

(ref . 8) 
N=4, M=16 

Lan 

(ref. 2'i 
N=8, M=30 

VLM 

(ref. 1) 
N=6, M=40 

Wagner 
from 
ref. 2 


Overall Values 

^L 

2.4732 

2.4744 

2.4707 

2.5239 

2.4778 

a 






-^M 

a 

0.5187 

0.5182 

0.5173 

0.5334 

0.5180 

X /C 
ac' 

0.2097 

0.2094 

0.2094 

0.2113 

0.2091 

1/e 

1.0007 

1.0007 

1 . 0050 

1.0018 

1.0005 


0.9951 

1.0108 

1.0022 

0.9764 

1.0172 


Values of Co/C, Values of x„ /c 

X, L ac 


2y/b 

Present 

NLR 

2y/b 

Present 

NLR 

0 

1.2543 

1.2543 


0 

0.2200 

0.2199 

0.1951 

1.2331 

1.2331 


0.1951 

0.2187 

0.2187 

0.3827 

1.1692 

1.1692 


0.3827 

0.2150 

0.2149 

0.5556 

1.0625 

1.0625 


0.5556 

0.2087 

0.2085 

0. .071 

0.9137 

0.9137 


0.7071 

0.1999 

0.1996 

0.8315 

0.7257 

0.7257 


0.8315 

0.1896 

0.1886 

0.9239 

0.5045 

0.5044 


0.9239 

0.1798 

0.1773 

0.9808 

0.2588 

0.2587 


0.9808 

0.1731 

0.1685 
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x/c = (1 - cos 0)/2 

0 Vortex Position: 0j^ = (k = 1 

X Control Points : 0^ = (i = 1, 

Figure 1.- Lan's vortex arrangement for airfoils 




Spanwise Control Points 


Rectangular Wing for M 


Figure 3.- Arrangement of trailing vortices for Prandtl’s 
lifting-line equation. 
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Figure 5.- Effect of number of trailing vortices on induced 
drag from Prandtl's lifting-line theory, A = 2tt. 
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Figure 6.- Arrangement of vortex lattice for lifting-surface equation 
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Figure 8.- Effect of vortex-lattice arrangement for 
rectangular planfonns. 


